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ABSTRACT 



Context. B[e] supergiants are known to possess geometrically thick dusty disks. Disk-forming wind models in the literature have, however, 
been found to be insufficient in reproducing the observed dust emission. This problem arises due to the severe assumption that, as for classical 
Be stars, the near-infrared excess emission originates in the disk. Modeling of the free-free and free-bound emission therefore results in an 
upper limit for the disk mass loss rate as well as for the disk opacity. Dust condensation in the disk can thus severely be hampered. 
y—^ Aims. In order to overcome the dust formation problem, and based on our high-resolution optical spectroscopy and model results, we propose 
\^ a revised scenario for the non-spherical winds of B[e] supergiants: a normal B-type line-driven polar wind and an outflowing disk-forming 
wind that is neutral in hydrogen at, or very close to the stellar surface. 

Methods. We concentrate on the pole-on seen LMC B[e] supergiant R 126 and calculate the line luminosities of the optical [Oi] emission 
lines and their emergent line profiles with an outflowing disk scenario. In addition, we compute the free-free and free-bound emission from a 
"j '' line-driven polar wind and model the spectral energy distribution in the optical and near-infrared. 
O Results. Good fits to the [Oi] line luminosities are achieved for an outflowing disk that is neutral in hydrogen right from the stellar surface. 

Neutral thereby means that hydrogen is ionized by less than 0.1%. Consequently, the free-free and free-bound emission cannot (dominantly) 
^ arise from the disk and cannot limit the disk mass loss rate. The hydrogen neutral outflowing disk scenario therefore provides an ideal 
. . environment for efficient dust formation. The spectral energy distribution in the optical and near-infrared range can be well fitted with the 
^ stellar continuum plus free-free and free-bound emission from the polar line-driven wind. Our modeling further delivers minimum values for 
A^disk S 2.5 X lO"' Mgyr"' and for the density contrast between equatorial and polar wind of ~ 10. 
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1. Introduction 

The spectra of B[e] supergiants show the so-called hybrid char- 
acter which is defined by the co-existence of a line-driven po- 
lar wind and a high density but low velocity equatorial wind 
(Zickgraf et al. ll985l) . The latter is assumed to form a disk-like 
structure. This disk is the location of the low-ionized metals, of 
the sometimes observed CO and TiO molecular emission bands 
(McGregor et al. IT9881 [19891 Zickgraf et al. [19891 ) as well as 
of the hot dust, pronounced in the mid-IR excess emission (e.g. 
Zickgraf , 1992) . The disk formation mechanism, as well as the 
disk structure, are however still unsolved problems. 

Porter ( I2003I I has studied the possible nature of these disks 
by modeling their spectral energy distribution (SED). He found 
that neither an outflowing disk-forming wind nor a Keplerian 
viscous disk in their simplest form can easily account for the 
observed free-free and dust emission self-consistently. 

Send offprint requests to: M. Kraus, 
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The optical spectra of B [e] supergiants exhibit strong [Oi] 
emission (see e.g. Kraus & Borges Fern andes 120051 Kraus et 
al. 120061 and Fig.[T]i. Besides other typical features like strong 
Balmer emission. Fen and [Feii] emission, and mid-IR excess 
emission, the [Oi] emission lines are one of the main character- 
istics of stars with the B[e] phenomenon (Lamers et al. |1998l l. 
The fact that Oi and Hi have about equal ionization potentials 
requests that the [Oi] emission region must be neutral in hydro- 
gen. Its most plausible location is therefore the high-density 
disk. 

Test calculations for an outflowing disk scenario, that were 
aimed to reproduce the observed strong [Oi] line luminosities, 
have been performed for some Magellanic Cloud B[e] super- 
giants (Kraus & Borges Fernandes 12005! Kraus et al. I2006I I. 
These computations revealed that, in order to keep the disk 
mass loss rates at reasonable values, the disk material must be 
neutral in hydrogen either at or at least very close to the stellar 
surface, resulting in a hydrogen neutral disk. 
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Fig. 1. The [Oi] lines in the FEROS spectrum of R 126. They 
all show a redshift of about 260 km s"' caused by the radial 
velocity of the LMC. 



In this paper we present model results for the Large 
Magellanic Cloud (LMC) B[e] supergiant R 126. We use a sim- 
ple outflowing disk scenario and request that the disk is neutral 
in hydrogen right from the stellar surface. We show that with 
such a model we can explain observations like the strong [Oi] 
line luminosities. In addition, this scenario provides a way out 
of the dust formation problem mentioned by Porter (I2003I I. The 
paper is structured as follows: in Sect.|2]we give details on the 
test star and on our observations. Sect.[3]introduces the problem 
of the nature of the B[e] supergiants disks. We emphasize the 
need of a revised disk model and propose the hydrogen neutral 
disk as a reasonable scenario. In Sect.|4]we present our hydro- 
gen neutral disk model and calculate the emerging [Oi] Une 
luminosities as well as their line profiles. In Sect. |5] we fit the 
SED in the optical and near-IR range by calculating the free- 
free and free-bound emission from the polar wind. In addition, 
we determine the density ratio between equatorial and polar 
wind components and show that our model results in a lower 
limit. The validity of our assumptions is discussed in Sect.|6] 
and the conclusions are summarized in Sect.]?] 

2. The LMC B[e] supergiant R 126 

The LMC B[e] supergiant R 126 (HD 37974, LHA 120-S 127) 
is supposed to have a pole-on seen disk. This orientation has 
the following advantages (see Fig.|2]i: (i) it guarantees that we 
see all of the [Oi] emission and none is hidden by the star or ab- 
sorbed by the dusty disk because the [Oi] emission originates 
from regions close to the star where no dust has yet formed, 
and (ii) it simplifies the analysis due to the high degree of sym- 
metry. R 126 has also been studied by Porter (I2003I I so that we 
can directly compare our results. The stellar parameters of this 
object have been derived by Zickgraf et al. (119851) who found: 
Teff = 22 500 K, U = 1.2 X 10'' Lq, = 72 Rq, and an inter- 
stellar extinction value of E{B -V)- 0.25. These values will 
be used for our calculations throughout the paper. 

We obtained high- and low-resolution optical spectra of 
R 126 at the ESO 1 .52-m telescope (agreement ESO/ON-MCT) 
in La Silla (Chile) using the Fiber-fed Extended Range Optical 
Spectrograph (FEROS) and the Boiler & Chivens spectrograph. 



respectively. FEROS is a bench-mounted Echelle spectrograph 
with fibers, that cover a sky area of 2 arcsec, located at the 
Cassegrain focus with a wavelength coverage from 3600 A to 
9200 A and a spectral resolution of R = 55 000 (in the region 
around 6000 A). It has a complete automatic on-line reduction, 
which we adopted. The spectra were taken on December 18, 
1999, with an exposure time of 4500 seconds. The S/N ratio in 
the 5500 A region is approximately 100. 

The low-resolution Boiler & Chivens (B&C) spectrum was 
taken on October 31, 1999, with an exposure time of 900 sec- 
onds and a slit width of 4 arcsec. The instrumental setup em- 
ployed provides a resolution of ~ 4.6 A in the range of 3800- 
8700 A. In the 5500 A continuum region, the S/N ratio is aprox- 
imately 200. The B&C spectrum was reduced using standard 
IRAF tasks, such as bias subtraction, flat-field normalization, 
and wavelength calibration. We have done absolute flux cali- 
bration using spectrophotometric standards from Hamuy et al. 
(fT994l) . 

Equivalent widths and line intensities in both linearized 
spectra have been measured using the IRAF task that computes 
the line area above the adopted continuum. Uncertainties in our 
measurements come mainly from the position of the underly- 
ing continuum and we estimate the errors to be about 20 %. The 
emerging [O i] lines resolved within the FEROS spectrum are 
shown in Fig.[T] 

3. On the nature of B[e] supergiant stars' disks 

The nature of the B[e] supergiant stars' disks is a long-standing 
problem. Recently, Porter (i2003i) investigated the possibility of 
dust formation in the disk of the B[e] supergiant R 126 for 
two different approaches: an outflowing disk-forming wind, 
and a Keplerian viscous disk. He found that both models failed 
in reproducing the observed dust and free-free emission self- 
consistently. And he suggested a way out of this by allowing 
for substantial alteration especially of the exponent of the radial 
density distribution from p ~ r"^ to a considerably flatter one 
of r"^-^. The advantage of a flatter density profile is that at ev- 
ery location in the disk a higher density (and hence opacity) is 
maintained, allowing for more efficient dust condensation and 
therefore for an enhanced dust emission over the disk. Such a 
modification finally resulted in a good fit to the SED of R 126 
(see Fig. 6 of Porter l2003l . but for the price that the disk density 
distribution now no longer corresponds to a constant velocity 
outflow (which requests p ~ r^^). A possible interpretation of 
such a flatter density profile might be that the disk slows down 
with distance, a scenario that lacks any observational evidence. 

To test Porter's modified outflowing disk scenario, we took 
the parameters of his best-fit model and calculated the emerg- 
ing [Oi] line luminosities. These turned out to be at least ?l fac- 
tor of 50 lower than the values we observed. This means that 
we need a much higher density within the [Oi] line forming 
region than can be provided by Porter's disk model. 

The most severe limitation in his model calculations was 
the assumption that, as in the case of classical Be stars, the free- 
free emission arises in the high density disk, while contribu- 
tions from the polar wind are negligible. Fitting the near-IR part 
of the SED with free-free emission from the outflowing disk 
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Fig. 2. Sketch of the pole-on seen disk around R 126 (not to 
scale). The disk (dark grey shaded) is assumed to be neutral 
in hydrogen right from the stellar surface, giving raise to [Oi] 
emission. Further out, molecules and dust will form. The disk 
material moves outwards with Vout, resulting in an observable 
line-of-sight velocity (v'los). Free-free and free-bound emission 
is restricted to the (polar) wind (light grey shaded). The sys- 
tem is axially symmetric (around z-axis) and symmetric with 
respect to the disk mid-plane. 



therefore determines the disk mass loss rate. Consequently, this 
value is an upper limit for the disk density and hence opacity. 
Dust can only condensate at distances with reasonable values 
of temperature and opacity. An upper limit for the disk den- 
sity, defined by the free-free emission, can therefore severely 
hamper efficient dust formation in the disk, which turned out 
to be the major reason why Porter (I2003I I could not fit the ob- 
served dust emission of R 126 with the original outflowing disk 
scenario. 

The detection of strong [Oi] line emission in the spectra 
of e.g. B[e] supergiants requests that there must be a rather 
high density disk region which is neutral in hydrogen {T % 
lOOOOK), but still hot enough (6000 < T [K] < 8000) for ef- 
fective excitation of the levels in Oi. Fkst test calculations that 
were aimed to reproduce the [Oi] line luminosities with an out- 
flowing disk scenario have been performed by Kraus & Borges 
Fernandes (2005) and Kraus et al. (I2006I I. These computations 
showed that, in order to keep the disk mass loss rates at reason- 
able values, the disks must be neutral in hydrogen either at, or 
at least very close to the stellar surface. 

The existence of hydrogen neutral disks around B [e] super- 
giants is additionally supported by recent ionization structure 
calculations in non-spherically symmetric winds. These com- 
putations have shown that the wind material in the equatorial 



plane can indeed recombine at, or close to the stellar surface 
if either the equatorial mass flux is enhanced compared to the 
polar one (Kraus & Lamers l2003l l. or even decreased in com- 
bination with a reduction in surface temperature due to rapid 
rotation of the star (Kraus 120061 Kraus & Lamers I2006I I. 

Having a hydrogen neutral disk around a luminous object 
is therefore not as unrealistic as it might seem at first glance. 
In fact, it has several striking advantages: if the disk is pre- 
dominantly neutral, the free-free and free-bound emission can 
no longer be generated in the disk, but must (mainly) originate 
within the polar wind. The disk mass loss rate is therefore not 
strictly linked to the free-free emission causing the near-IR ex- 
cess. Hence, the disk mass loss rate might be (much) higher 
than the upper limit used by Porter (I2003I I. A higher disk den- 
sity is certainly needed to reproduce the observed [Oi] line lu- 
minosities. And finally, a higher disk density will increase the 
disk opacity and will allow for more efficient dust condensation 
in the disk, needed to reproduce the observed IR excess. 

In the following we will investigate and test the suggested 
model scenario of a hydrogen neutral disk around a luminous 
object by analysing in detail the observed [Oi] lines from the 
B[e] supergiant R 126. 

4. The (hydrogen) neutral disk 

4.1. The disk geometry 

The disk is assumed to be wedge-shaped (see sketch in Fig.|2|l 
and we use the same value of a = 10° for the constant disk 
half opening angle as Porter (120031 ). For simplicity we assume 
that the disk has a constant outflow velocity within the [Oi] 
emitting region. The influence of a velocity distribution on the 
final results is discussed in Sect.|6] 

In a disk-forming wind of constant outflow velocity, the hy- 
drogen particle density distribution at all latitudes within the 
disk can be written in the form 



tinir) = nH(^in)^ ■ 

The parameter nH(^in), given by 

\ I j RA^ ^'m.disk 
nH(^in) = — 



jUOTh \R 



(1) 



(2) 



with the disk mass flux, Fm.disk, and the (constant) disk out- 
flow velocity, Vout (see Fig.|2]l, describes the density at the inner 
radius, Rin, of the [Oi] emitting region. This inner radius corre- 
sponds to the hydrogen recombination distance in the disk and 
is in principle a free parameter in our calculations, but here we 
use Rin - Rt since this defines a lower limit to the disk den- 
sity needed to reproduce the [Oi] line luminosities. The valid- 
ity of our assumption that the disk is neutral in hydrogen right 
from the stellar surface is discussed in Sect.|6] The oxygen and 
electron density distributions are parameterized in terms of the 
hydrogen density, i.e. 



no(r) = qon}i(r) and n^^ir) = qenii(r) 



(3) 



The parameters qo and q^ describe the relative abundances of 
oxygen and electrons with respect to hydrogen. For qo we use 
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Fig. 3. Sketch of the five lowest energy levels within Oi from 
which the forbidden lines arise. The transitions of our interest 
are indicated. 



a typical LMC abundance value of 1/3 solar, with a solar oxy- 
gen abundance of ^o.soiai - 6.76 x 10"^ (Grevesse & Sauval 
I1998I I. Since we do not know the ionization fraction of the disk 
a priori, we keep as a free parameter. 

4.2. The [Oi] lines and the disk ionization 

The forbidden emission lines arise from the five lowest en- 
ergy levels of the Oi atom that are excited via collisions with 
free electrons. These levels are shown schematically in Fig. [3] 
For simplicity, we labeled them from 1 to 5 corresponding to 
their increasing energy. The three lines resolved in our high- 
resolution optical spectra are indicated; they belong to the tran- 
sitions 5^ 4 (5577 A), 4 1 (6300 A), and 4 ^ 2 
(6363 A). 

While the lines 6300 A and 6363 A arise from the same 
upper level, which means that their line ratio is independent 
of the density and temperature but is determined purely by 
atomic parameters, the situation is different for line ratios with 
the 5577 A line. This line arises from the highest energy level 
in our 5 -level atom (Fig.|3]l. It is therefore obvious that the 
5577 A line will only become strong if level 5 is populated effi- 
ciently, either due to a huge number of available free electrons 
or by electrons with rather high energy. Line ratios with the 
5577 A line are thus sensitive tracers for the temperature and 
ionization fraction in the disk. In our spectra, the 6300 A line 
is the one with highest luminosity. We will therefore use the 
6300 A/5577 A line luminosity ratio to investigate density and 
temperature in the disk around R 126. 

Forbidden lines are optically thin. The emissivity at any 
location r in the disk of a transition n — > m is given by 
jn,mif) - /!V'n,mnn('"Mn,m where «„('") IS the level population 
of the upper level, n, at this location. The level populations are 
calculated by solving the statistical equilibrium equations in 
our 5-level atom. Collision parameters are taken from Mendoza 
( I1983I I and atomic parameters from Wiese et al. (.1966 ) and 
from Kafatos & Lynch ( |19801 l. The line luminosity finally fol- 



lows from integration of the emissivity over the emitting disk 
volume. 

Due to the nearly equal ionization potentials of O and H, 
the [Oi] emission must arise in those parts of the disk in which 
hydrogen is predominantly neutral. This does not mean that hy- 
drogen must be 100% neutral, since effects like collisional ion- 
ization that might take place in the high density disk could keep 
the hydrogen material ionized at a (very) low level, providing 
free electrons to collisionally excite Oi. Further sources of free 
electrons in the disk are elements like Mg, Si, Fe, Al, etc. with 
ionization potentials (much) lower than 13.6eV. These metals 
will remain ionized at a certain fraction even in a completely 
hydrogen neutral environment. The ionization fraction, q^, in 
the disk is therefore the sum of all free electrons provided by 
hydrogen and the metals, i.e. - qu+ + q'Metais+- 

To get a handle on the possible individual contributions to 
the total number of free electrons in the disk, we make a first 
guess. If all possible elements with ionization potential lower 
than 13.6eV would be fully singly ionized (including also C 
with^ = 11.26eV andCl with^ = 12.97 eV), then a maximum 
number of free electrons of ^Metais+ ^ 1.57 x lO"'* would be 
provided by the metals in the disk. This number accounts for 
an LMC metallicity of 1/3 solar On the other hand, if only 1 % 
of hydrogen in the disk remains ionized, the ionization fraction 
would be q'H+ - 10"^ which is at least a factor 100 higher than 
what can be expected from the metals. Our assumption of a 
hydrogen neutral disk around R 126 therefore requests that the 
total disk ionization fraction should be q^ < 10"^. 

To derive the ionization fraction in the disk of R 126 we 
calculated the line luminosities and the 6300 A/5577 A line 
ratio for a large range of values and for different temper- 
atures. These calculations were done in several steps: for a 
given ionization fraction and temperature we first varied the 
density parameter, F^ disk/vout (see Eq.|2]l, in order to repro- 
duce the observed line luminosity of the 6300 A line. This 
density parameter is plotted as a function of in the lower 
panel of Fig.|4] With the proper density parameter, we then cal- 
culated the line luminosities of all [Oi] lines and derived the 
6300 A/5577 A line ratio. This line ratio is plotted as a function 
of qs in the upper panel of Fig.|4] In this plot we also included 
the observed line ratio, shown as the dotted line. 

Our calculations cover a temperature range between 7500 K 
and 9000 K. The curves in the upper panel of Fig. |4] all show 
the same trend of a decrease in line ratio with increasing ion- 
ization fraction. This effect is clear, because the number of free 
electrons determines the collisional excitation rate of the levels. 
The more electrons, the higher the probability for collisional 
excitation even of the level with highest energy (i.e. level 5) 
from which the 5577 A line arises. The second obvious effect 
is the decrease in line ratio with increasing temperature, since 
the temperature determines the kinetic energy of the free elec- 
trons. For low temperatures, the energy of the free electrons 
is not high enough to effectively excite level 5, reducing the 
5577 A line emission and hence increasing the line ratio. 

For the three temperature values used to calculate the theo- 
retical curves in Fig.|4]we find agreement with the observed line 
ratio for the following parameter combinations: for 9000 K we 
find an ionization fraction of 10"^, and for 8000 K we find an 
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Fig. 4. Upper panel: Variation of the line ratio with disk ioniza- 
tion fraction. The three curves are for the different temperatures 
as indicated. The dotted line represents the observed line ratio 
from R 126. It indicates a disk temperature in the range 9000 K 
to about 8000 K with corresponding ionization fractions of 
10"^ and 4 x 10""*, respectively. Lower panel: Density parame- 
ter of the disk found from fitting the observed 6300 A line lu- 
minosity over a huge range of ionization fractions. For details 
see text. 

ionization fraction of 4 x 1 For a temperature of 7500 K, we 
cannot reproduce the observed line ratio. 

Our calculations show that when accounting for the line 
luminosity and line ratio observed from R 126, the ionization 
fraction is rather low for high temperatures, while an ion- 
ization fraction of up to 10% might exist for a temperature 
as low as 7800 K. This behaviour is only strictly valid for a 
given observed line luminosity. In order to restrict the range 
of valid or reasonable disk parameters, we take account of the 
fact that the ionization fraction can only increase with increas- 
ing temperature, since collisional ionization will act more ef- 
ficiently. This indicates that for the disk around R 126 we can 
exclude high temperatures 9000 K) as well as low temper- 
atures 8000 K), and we conclude that the [Oi] line form- 
ing region in the disk of R 126 has a temperature range of 
~ 8000 . . . 8500 K with a corresponding ionization fraction 
range of 5 x 10 ^ ~ 5 x 10"^. This range in ionization 
fraction states that less than 0.05% of hydrogen in the disk is 
ionized. We can therefore conclude that the disk around R 126 
is indeed neutral in hydrogen. 

Fig.|5] shows the increase in [Oi] line luminosities calcu- 
lated for a disk of 8000 K and an ionization fraction of 4 x 10""*. 
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Fig. 5. Calculated increase in [Oi] line luminosities from the 
hydrogen neutral disk around R 126. The observed, extinction 
corrected values are indicated with arrows. The observations 
have been taken with an aperture radius of 2", corresponding 
to an outer edge of the observed wind of about 1.5 x 10^ R© 
at the distance of the LMC. This value is used as the extent 
of the line luminosity integration, but the [Oi] line luminosities 
clearly saturate already within about 500 R,. 

The lines saturate within a few hundred stellar radii. The ob- 
served values are indicated by the arrows to the right. 

4.3. Disk outflow velocity and mass loss rate 

For a temperature of 8000 K as derived in the previous sec- 
tion, we end up with a density parameter of Fm.disk/i'out - 
2.2 X 10"" gcm"^ for the disk of R 126 (see Fig.|4|i. To com- 
pare our results with those of Porter (120031 1. we need to convert 
this density parameter into a mass loss rate. For this, the out- 
flow velocity, Vout, of the disk is needed. For a pole-on seen 
disk with radially outflowing material, the maximum observ- 
able line-of-sight velocity, vios, is linked to the outflow velocity 
via v'los = Vout sin a (see Fig.|2]i. The range in possible vios val- 
ues over the different disk outflow directions above and below 
the disk mid-plane leads to symmetrical doppler shifts of the 
line center, while thermal and turbulent velocities as well as the 
spectral resolution of the instrument result in a Gaussian pro- 
file of the line. Thus, optically thin lines formed in a pole-on 
seen outflowing disk should be symmetric with respect to their 
laboratory wavelength. The observed [Oi] lines are indeed sym- 
metric (Fig.[T]l and their profiles show a Gaussian shape with a 
FWHMof ~ 9.2kms-i. 

To derive a possible outflow contribution from the observed 
lines, we first calculate the following test profiles: (i) a pure 
Gaussian line profile and (ii) a pure outflowing disk profile. 
The results are shown in the left panel of Fig.|6] The calcu- 
lated lines are normalized to their maximum value, and Vgauss 
and Vout are chosen such that they result in a FWHM value of 
9.2km s '. The pure outflowing disk (neglecting the influence 
of FEROS' spectral resolution) results in a flat-topped profile 
with steeply rising wings. The narrow observed [Oi] lines show 
a better agreement with a Gaussian line shape, and the fact that 
the observed [Oi] lines have widths that are only marginally 
broader than the velocity resolution of FEROS indicates that 
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Fig. 6. Normalized profiles of the [Oi] 6363 A line. Left: 
Theoretical profiles for either a pure Gaussian or a pure out- 
flowing disk profile. Middle: Fit of a Gaussian profile with 
FWHM of 9.2 km s ' resulting from FEROS' spectral reso- 
lution plus thermal motion of the gas. Right: Fit of the same 
Gaussian including an outflow component causing vios of either 
2 km s"' or 4 km s"' . While the first one still gives a reasonable 
fit, the latter one clearly results in a too broad line profile. 

the intrinsic FWHM of the lines can only be on the order of 2- 
3kms ^ Such a low intrinsic FWHM can be ascribed to pure 
thermal broadening. 

Since spectra taken with FEROS are not flux-calibrated, the 
reproduction of the line profiles can only give qualitative re- 
sults. We took exemplarily the line profile of the 6363 A line 
obtained with FEROS, corrected it for the LMC radial velocity 
component of 260km s ', and normalized the line to its maxi- 
mum intensity. This normalized line profile is what we intend 
to fit. 

We first used a pure Gaussian line profile with FWHM 
of 9.2 km s ' (resulting from FEROS' spectral resolution plus 
thermal broadening of the line) which results in a good fit (mid 
panel of Fig.|6l). 

Next we combined this Gaussian profile with an outflowing 
disk profile for different values of vios. This was done to see for 
which value of vios the line profile starts to alter. We found that 
a possible vios contribution must be smaller than ~ 2kms"' 
(right panel of Fig.|6l) resulting in an outflow velocity of Vout ~ 
ll.Skms ' for a - 10°. The value of a should be taken as 
an upper limit since the range of possible B[e] supergiant disk 
values has been estimated by Zickgraf ( 119921 ) to be a - 5-10°. 
A thinner disk will result in a higher outflow velocity for the 
same value of v'los (see Fig.|7ll. Nevertheless, even considering 
a thinner disk, the outflow velocity is considerably smaller than 
the values of 60-80 km s ' usually ascribed to the disks of B[e] 
supergiants. In the case of R 126, no good estimates for the real 
disk outflow velocity exist and only a detailed investigation of 
the line profiles from emission lines that might be connected 
with the disk will help to clarify this point. 

An outflow velocity of ll.Skms"' as derived above to- 
gether with the disk density parameter derived from the [Oi] 
line luminosity calculations and a disk opening angle of 20° re- 
sults in a total disk mass loss rate of Mdisk - 2.5 x 10"^ Moyr"' . 
This is about a factor of 10 higher than the value used by 




5 ^10 

disk hal" opening angle a [ ] 

Fig. 7. Disk outflow velocity as function of disk half opening 
angle for vios = 2kms"' as derived from the [Oi] line profile 
fitting. 

PorteiQ(|2003]l, and well in agreement with the postulated need 
for a higher disk density in order to account for the observed 
[Oi] line luminosities of R 126. 

Based on recent observations with the Spitzer Space 
Telescope Infrared Spectrograph, Kastner et al. (I2006t derived 
a total dust mass of the disk around R 126 of ~ 3 x 10"^ M©, 
and they claimed that the dust is contained between 120AU 
and 2500 AU. In order to compare this dust mass with our de- 
rived disk mass loss rate, we assume a gas to dust ratio of about 
100 and convert the total mass into a disk mass loss rate for our 
outflowing disk model. We end up with a density parameter 
of /^m.disk/i'out - 3 X lO^'^gcm""' and a disk mass loss rate of 
A^fdisk - 3.4x lO^'^Moyr"'. This is about 170 times higher than 
the value used by Porter (2003 ) and about 13 times higher than 
our value, and confirms that the disk around R 126 must indeed 
be more massive than previously thought. 

In addition, the inner edge of the dusty disk of 120 AU, 
which corresponds to about 360 R,, lies within our [Oi] satura- 
tion region. This confirms that our value for the disk mass loss 
rate is indeed a lower limit. 

5. The polar wind 

With our assumption of a hydrogen neutral disk right from 
the stellar surface, the free-free and free-bound emission is re- 
stricted to the ionized wind parts (see Fig.|2]i. At all latitudes 
above and below the equatorial disk we assume that the star 
has a normal B-type line-driven wind, which we refer to as the 
polar wind. Its radial density distribution, following from the 
equation of mass continuity, is 

1 \^ /^m.pol 

nnir) = — — — , (4) 

pnin \ r I v(r) 

' Note that the parameter used by Porter does not describe the 
total disk mass loss rate. Instead, his value of Md = 10"'' Moyr"' 
needs to be multiplied by the solid angle of the disk and divided by 
An. The total disk mass loss rate in his models is therefore Mask - 
2x lO-^Moyr-i. 
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where F^ ^oi is the mass flux of the polar wind. The velocity 
increase in line-driven winds can be approximated with a /3- 
law with /3 typically in the range of 0.8-1.0 for hot star winds 
(see e.g. Lamers & Cassinelli |1999l ). We assume that the wind 
temperature is constant which allows for a simplified treatment 
of the intensity calculation. At every point along the impact 
parameter, ( (see Fig.|2]l, the intensity is given by 

/, = B,(r)(i-e-(^>) . (5) 

The optical depth Ty(^) = J Ky(^, s)ds is defined by the line-of- 
sight integral (with the line-of-sight perpendicular to the disk 
mid-plane) over the absorption coefficient, Ky, defined by 

Kyi^, s) - 3.692 X 10*^ ^^^^ (l - e-^) (gg + gf,) , (6) 

and integration limits depending on the geometry of the system 
(Fig.|2]i. The parameters gg and gfb are the gaunt factors for 
free-free and free-bound processes, respectively, and have been 
calculated (e.g. by Kraus '2000) over a large frequency range 
and for different temperatures. Due to the symmetry and the 
assumed constant wind temperature, the total observable flux 
density of the free-free and free-bound emission is given by 

Fy^^ / 5v(r)(i-e-<f>)^fl'^ (7) 



where ^^ax is the outer edge of the ionized wind and d is the 
distance to the object. Due to the pole-on view, the dusty disk 
around R 126 might cause absorption of the free-free and free- 
bound emission from the rear parts of the wind (i.e. below the 
disk). However, as stated by Kastner et al. (12006b . the dust is 
expected to exist only at large distances, i.e. beyond 120 AU. 
For the wind densities of our interest, we found that the free- 
free and free-bound emission is generated much closer to the 
star, and any absorbing influence of the dust can be neglected. 

To calculate the electron density at every point in the wind, 
we need to know the polar mass loss rate and terminal veloc- 
ity. Literature values for R 126 are still not well constrained 
and vary over a huge range, depending on the geometry and 
the method used by individuals to derive them. As an example 
serves the terminal velocity for which Bjorkman (1998) gives 
650 km s ' while Zickgraf et al. (fT985 i) found 1 800 km s ' . Due 
to this uncertainty by almost a factor of 3 we calculate the free- 
free and free-bound emission by varying the density param- 
eter, i.e. the ratio F^ po\/Voo, instead of the individual param- 
eters F^ poi and Vco- The total continuum, including the stel- 
lar emission, is then reddened with the extinction of R 126, 
and compared with photometric data from the literature. A 
good fit to the SED in the optical and near-IR is achieved for 
f m,poi/voo ^ 2.0 X 10"'^ g cm-^^ and /3 = 0.8 (see Fig.®. 

The density parameters for the equatorial and polar wind 
as summarized in Table[T] result in a the density contrast on 
the order of Peq/Ppoie - 10. This value turns out to be a lower 
limit, because all parameters are chosen such that any change 
(like (i) a smaller disk opening angle, (ii) the existence of a 
small ionized inner disk part, (iii) a vertical density distribution 
within the disk due to a latitude dependence of the parameter 
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Fig. 8. Calculated SED of R 126 consisting of the stellar contin- 
uum and the free-free and free-bound continuum from the polar 
wind. The theoretical spectrum is reddened to the R 126 value. 
Photometric data are from Zickgraf et al. ( 119851 open circles), 
and from the public databases DENIS (squares), 2MASS (tri- 
angles), and MSX (filled circles). Not included here is the mod- 
eling of the dust emission to account for the mid-IR excess. 

Table 1. Model parameters for the polar and equatorial (i.e. the 
disk-forming) winds. Please note that the disk parameters only 
hold for the [Oi] Une forming region. 



Wind 






9e 




[K] 


[gem"'] 




Polar 


10 000 


2.0 X 10"'- 


1.0 


Disk-forming 


8 000 


2.2 X 10"" 


4.0 X lO"'' 



nH(^in), (iv) a cut-ofF in [Oi] line luminosities before saturation 
is reached) will immediately request a higher disk density in 
order to reproduce the observed [Oi] line luminosities. 

With a disk density parameter of 10 times the polar one, but 
a disk ionization fraction of less than 10"^, the radial electron 
density distribution within the disk is at least a factor of 100 
lower than the one in the polar wind. A contribution of the disk 
to the free-free and free-bound continuum is therefore indeed 
negligible, justifying our assumption that the free-free and free- 
bound emission is generated purely within the polar wind. 

6. Discussion 

For our description of the disk model in Sect. l4.1l we made two 
severe assumptions, namely we requested that (i) the disk is 
neutral in hydrogen right from the stellar surface, i.e. - R», 
and that (ii) the disk has a constant outflow velocity. We will 
now discuss the validity of these assumptions. 

A hydrogen neutral disk right from the stellar surface. 
We requested that the disk around the rather hot central star 
is predominantly neutral in hydrogen right from the stellar 
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surface. Recent ionization structure calculations by Kraus & 
Lamers (I2003I I in non-spherical winds of B[e] supergiants have 
shown that in the equatorial region hydrogen can indeed re- 
combine close to the stellar surface, resulting in a geometri- 
cally thick hydrogen neutral disk. To quantify this, we refer to 
their model C which was calculated for a star with the same ef- 
fective temperature as R 126. Further parameters given in their 
Table 1 were an electron temperature of lOOOOK, and an equa- 
torial surface density of 3.3 x 10'^ cm"-'. With these parameters, 
Kraus & Lamers (2003 ) found that hydrogen recombines in the 
equatorial plane at a distance of r ^ 1.0004R,. 

Our model calculations for the [Oi] line luminosities re- 
vealed a disk density parameter of about 2.2 x 10"" gem"-' 
(see Table[TJ which results in an equatorial density on the stel- 
lar surface of riuiR*) - 10'-' cm"^. This density is about 3 
times higher than the one used by Kraus & Lamers (2003). 
Consequently, recombination of the disk material will happen 
even closer to the stellar surface. It is therefore reasonable to 
adopt that the disk around the B[e] supergiant R 126 is neutral 
in hydrogen right from the stellar surface. 

The constant disk outflow velocity. In order to describe 
the disk density distribution with only one free parameter, i.e. 
the density parameter nH(^in) ~ ^'m.disk/i'out - const, we set the 
disk outflow velocity constant within the [Oi] line forming re- 
gion. Fig.|5]shows that the line luminosities of the [Oi] lines sat- 
urate within about 500/?*; the 5577 A line saturates even within 
50 Rf Let us assume that the disk-forming wind also has a ve- 
locity distribution according to the y6-law. Such a velocity law 
then mainly influences the formation of the 5577 A line. This 
line is created mainly in those wind parts in which the wind is 
still accelerating, while the other two [Oi] lines are unaffected 
from a y6-type velocity law. Their luminosities are produced at 
distances at which the velocity has reached already its terminal 
value. 

According to the equation of mass continuity, Eq. (|4]i, a 
lower velocity at a certain distance results in a higher wind den- 
sity and therefore in a higher level population. Consequently, 
the luminosity of the 5577 A line arising in the accelerated 
wind region is somewhat enhanced compared to the lumi- 
nosity created in an outflow of constant (i.e. terminal) ve- 
locity. A higher line luminosity, however, results in a lower 
6300 A/5577 A line ratio. Therefore, in order to fit the observa- 
tions, a slightly different combination of disk ionization frac- 
tion and density parameter will be necessary. 

We calculated the disk ionization fraction and density pa- 
rameter for an outflowing disk with a temperature of 8000 K 
and a velocity law with /3 - 0.8. From our fitting procedure 
we found the following results: ^ 10"^ and Fm,disk/voo - 
1.7 X 10"" gem"-'. Even though the ionization fraction is a 
factor of 2.5 higher than in the case of the constant outflow 
velocity (see Table[T]i, the total amount of ionized hydrogen 
in the disk is still only 0.1%. In addition, the total disk mass 
loss rate found with the y6-type disk wind scenario is now 
^disk - 2 X 10"^ Moyr"'. This value is only slightly lower than 
the value derived with the constant outflow velocity scenario in 
Sect. l4.3l Therefore, the conclusion of a hydrogen neutral disk 
with a density ratio of about a factor of 10 between equatorial 
and polar wind remains valid. 



We want to emphasize that not much is known about the 
real velocity distribution of the outflowing disks around B[e] 
supergiants. Since a y6-type velocity law does not drastically 
alter the results, the assumption of a constant outflow velocity 
as a first guess seems therefore to be a reasonable approach. 

7. Conclusions 

This paper is not aimed to present a detailed analysis of R 126 
and to derive its wind and disk parameters with high accuracy, 
but to present a somewhat revised model for the outflowing 
disk scenario. We suggest that the disks around B[e] super- 
giants are neutral in hydrogen at, or close to the stellar surface. 
Indications for such a model come from observations as well 
as from theory and are here briefly summarized: 

- Former investigations of the disks around B[e] supergiants 
assumed that the observed near-IR excess is caused by free- 
free emission from the disk, as it is the case for classical Be 
stars. A contribution from the polar wind was thought to 
be neglibible. Fitting of the near-IR excess with free-free 
emission arising purely in the outflowing disk delivers a 
disk mass loss rate that defines an upper limit to the disk 
(dust) density. The resulting dust emission was however 
found to be too low to account for the observed IR excess 
(Porter "2003^. 

- Porter (2003 ) tried to solve this dust emission problem by 
allowing for a flatter disk density distribution. However, his 
modified disk model for the B[e] supergiant R 126 is not 
able to reproduce the observed [Oi] line luminosities, un- 
derestimating them by about a factor of 50! To account for 
the [Oi] line luminosities a (much) higher disk density is 
needed. 

- Due to the nearly equal ionization potentials of Oi and Hi, 
the existence of strong [Oi] emission further requests that 
in the line forming region hydrogen must be predominantly 
neutral. 

- Recent ionization structure calculations in non-spherically 
symmetric winds of B[e] supergiants have revealed that the 
wind material indeed recombines in the equatorial plane al- 
ready at, or very close to the stellar surface, resulting in a 
hydrogen neutral equatorial disk or zone (see e.g. Kraus & 
Lamers, |2003l|2006l Kraus|2006ll. 

The most important consequence of such a hydrogen neutral 
disk is the fact that the disk can no longer be considered as the 
formation location of the free-free emission causing the near- 
IR excess. This emission must then mainly be produced in the 
line-driven polar wind. Consequently, the fitting of the free- 
free emission does not result in an upper limit to the disk mass 
loss rate anymore. Instead, the disk mass loss rate, and hence 
the disk density, can be (much) higher, providing an ideal envi- 
ronment for the [Oi] line emission as well as for efficient dust 
condensation to account for the observed IR excess. 

We tested the scenario of a hydrogen neutral disk for the 
LMC B[e] supergiant R 126 by modeling the line luminosities 
and line ratios of the [Oi] emission lines resolved in our high- 
resolution optical spectra. The parameters derived for the disk 
and wind of R 126 are the following: 
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- We found that the [Oi] 6300 A/5577 A line ratio is very sen- 
sitive to the ionization fraction in the disk. From fitting the 
observed line ratio we can conclude that hydrogen in the 
disk is ionized by less than 0.1%. This confirms that the 
disk is indeed predominantly neutral in hydrogen. 

- The disk mass loss rate of Mdisk ~ 2.5 x lO"*" Moyr"' found 
from our fitting is about a factor of 10 higher than the value 
used by Porter ( I2003I I. This is in good agreement with the 
postulated need for a higher disk density in order to fit the 
[Oi] lines, and also in good agreement with the total dust 
mass of about 3 x lO^-'M© as derived by Kastner et al. 
(12006 V A disk with (much) higher density provides a much 
better environment for efficient dust condensation than the 
disk model used by Porter ( I2003I I. 

- The near-IR excess is fitted with free-free and free-bound 
emission from the B-type line-driven polar wind. The re- 
sulting density contrast between equatorial and polar wind 
is on the order of 10. This value is found to be a lower limit. 

To summarize, based on a detailed investigation of the emerg- 
ing [Oi] emission lines we found that the disk around the B[e] 
supergiant R 126 is neutral in hydrogen right from the stellar 
surface. Since all B[e] supergiants show strong [Oi] line emis- 
sion, we postulate that the hydrogen neutral ouflowing disk sce- 
nario might also hold for the other members of the B[e] super- 
giant class. 
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